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I. INTRODUCTION 

Soon after achieving the Bose-Einstein condensation 
in ultracold dilute atomic gases [H, there came a surge 
of interest in degenerate quantum gas mixtures that has 
been continuing unabated up to the present days. Multi- 
component Bose condensates were initially realized 0, |3| 
by using two different hyperflne states of ^''Rb, later mis- 
cible and immiscible spinor condensates of ^"^Na were ob- 
served Q. Heteronuclear condensates of ^^K and *^Rb 
have been realized in a trap fsl] and in an optical lattice 
0. Very recently, different groups have used Feshbach 
resonances to engineer two-component Bose condensates 
with tunable inter- and intra-species interactions in *^Rb- 
^^K 0, s^Rb-'^iK 0, and 87Rb-85Rb [9] mixtures. 

The interest to multicomponent Bose systems is 
driven, particularly, by their rich physics: they provide 
an opportunity to study various phase transformations 
ranging from demixing of interpenetrating quantum liq- 
uids [lO, im to transitions between different magnetic 
states in spinor condensates 12, 13], and to a variety of 
yet unobserved quantum phase transitions predicted to 
appear in the presence of an optical lattice |14ll8| . Both 
in traps and optical lattices, the realization of highly 
anisotropic geometries ( "cigars" and "pancakes" ) is pos- 
sible [l9|, which provides an opportunity to study the 
above phenomena in lower-dimensional systems. 

The simplest phase transition in a Bose mixture is the 
phase separation (demixing) in a two-component system, 
which occurs if the inter-species repulsion overcomes the 
intra-species one. Existing theoretical studies of this in- 
stability are usually performed at the mean-field level 
by using coupled Gross-Pitaevskii equations [l3, [HI . It 
is well known that in lower dimensions mean-field ar- 
guments might become inapplicable in the low density 
regime In a one-dimensional Bose gas, with the 

decrease of the density, the healing length ^ becomes 
smaller than the average interparticle distance d, thus in- 
validating the mean-field approach. In two dimensions, 
the situation is more subtle: although the ratio ^/d does 
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not depend on the density and falls below unity only in 
the case of very tight two-dimensional confinement, in 
the dilute limit the effective coupling constant becomes 
strongly energy- and density-dependent [2l|, HI] . 

For one-dimensional (ID) multicomponent systems, a 
study beyond the mean-field approximation, based on the 
bosonization technique |23l]. is available [2J|. However, 
the approach of Ref. [2J] is constrained by the require- 
ment that the inter-species coupling is small compared 
to the characteristic bandwidth, which becomes too re- 
strictive for dilute systems and leads to a breakdown of 
bosonization at very low densities. 

The goal of the present work is to provide a stability 
analysis for dilute multicomponent Bose gases in low di- 
mensions, valid for any coupling strength. It is shown 
that strong renormalization of coupling constants, typi- 
cal for low-dimcnsional systems, can lead to a substantial 
change in the stability conditions compared to mean-field 
results, particularly in the case of inequivalent species 
(heteronuclear mixtures); the mean- field answers restore 
their validity with the increase of density. The outline of 
the paper is as follows: in Section |ll] I present the gen- 
eralization of the renormalization group (RG) approach 
to the multicomponent case, in Section IIIII I apply this 
general formalism to the simplest case of two bosonic 
species with density-density interaction, and Section ITVl 
illustrates the case of mutually convertible species on the 
example of spin-one bosons. Finally, Section fVl contains 
the discussion of the possible numerical and physical tests 
of the theory predictions. 



II. RENORMALIZATION GROUP APPROACH 
FOR A MULTICOMPONENT DILUTE BOSE GAS 

As a starting point for our analysis, we will choose the 
theory of multi-species bosonic field with pointlike two- 
body interactions, that effectively describes a system of 
atoms with typical momenta much smaller than the in- 
verse characteristic potential range (atom size). Effec- 
tive field theory in combination with the renormalization 
group analysis is a convenient tool j25l-[27| that is widely 
used [28 30] for the study of the cold atom systems. 

Consider the following continuum action describing N 
species of bosons at zero temperature in d spatial dimen- 
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sions, coupled by a general quartic contact interaction: 



An = / dr / (9r - Ma)V'a + 



2ma 



+ U 



al3a'j3' 



(1) 



Here ipa sue the fields describing bosonic particles with 
masses nia, a — 1, . . . N, and we have set h — 1. The 
interaction matrix for Bose fields satisfies the obvious 
symmetry conditions gafi,a'f}' = gpa,a'i3' = gai3,i3'a' = 
9a'fi' Q/3- The action ([T|) can describe a mixture in a con- 
tinuum as well as in an optical lattice [23 |. 

When all chemical potentials /Xq, vanish, the system 
is at a critical point. At this special point, the phys- 
ical picture is considerably simplified [3l|, |32|: there is 
no self-energy correction so that the full propagator just 
coincides with the free one Gaf3{k,uj) — Sap/{iio — e^""*), 
where e^"'' = k'^/2ma. This is easy to understand phys- 
ically, since when all /i^ = 0, the particle density is just 
zero. Interaction between the particles is, however, renor- 
malized due to multiple scatterings. The renormalized 
vertex F, defined at the fixed sum of external momenta 
Q, satisfies the following Bethe-Salpeter equation (illus- 
trated in Fig. [T]) 

^aP.-iSiQ) = gafl,fS — gal3,a'l3' fa'P'{Q)^a'l3' ,fS{Q), 
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The simple ladder form of the above equation is due to 
the fact that contributions from any diagrams containing 
closed loops vanish [31|J. 

Above the upper critical dimension d = 2, fap{Q) 
converges at Q — >■ and ([2]) simply yields the dressed 
(observable) interaction matrix. For d < 2 there is a 
singularity at Q — > 0: fapiQ) oc 1/Q for d = 1 and 
fap{Q) ln(Ao/(5) for d ~ 2, where Aq is the ultraviolet 
cutoff (for a system in an optical lattice, the parameter 
Aq has the physical meaning of a lattice cutoff, and in a 
continuum it has the sense of an inverse characteristic po- 
tential range) . Thus for d < 2 one should look at the RG 
fiow of the running coupling matrix r{l) = F((5)e*-2~''-'' 
with the change of the scale Q i— )■ Age"'. Defining the 
matrices 

Tal3,a'l3' — ap ,a' P' {mal3ma' {3')^^^ , 
Fafl.a'fl'iQ) — ^al3,a'l3'Rd{Qal3), (3) 

where ma/s = rriamf^ / {ma + w^g) are the reduced masses, 
Qap = 2Q{mampY/'^ / [ma + mp), and 
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one can write down the RG equation as 

,{r(i-Fr)-i} = 0. (4) 
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Further, using the identity AX = XA = A — 1 with 
A = (1 - A)-i and its derivative dX/dl = X{dA/dl)X, 
those equations can be rewritten in the following form: 



dV/dl + T{dF/dl)V = 0. 



(5) 



For further simplification, it is convenient to redefine the 
interaction matrix once again. In one dimension (d = 1), 
one can introduce F — F^^'^TF^^^; then, since F is a 
diagonal matrix that for — 1 is proportional to 1/Q, 
one obtains dF/ dl ^ F and Eq. ([5|) can be rewritten as 



dV/dl = F-r2 



(6) 



which is familiar from the one-component case [33j; the 
only difference is that the interaction is now a matrix. 

In two dimensions (c? = 2), it is convenient to define 
F = r/27r; then, using the fact that dF/dl = l/(27r) for 
d = 2, one can reduce ^ to 



dV/dl = -F2 



(7) 



Summarizing the above derivation, we see that the 
RG equations for the interaction matrix of the low- 
dimensional multicomponent Bose gas can be cast into 
the common form 



dT/dl = {2- d)T -T"^, 



(8) 



where 

ra/3,7i5(0 



^ap.,'ysil)y/ma0m^s/T^, d = 2. (9) 



Since the F matrices are symmetric, they can always be 
diagonalized by an appropriate orthogonal transforma- 
tion. It is worth noting that at the zero-density critical 
point the equations ([8]) are exact to all orders in the in- 
teraction F, similar to the one-component case (33j . 

The above derivation is nothing but a direct gener- 
alization of the well-known RG approach to the one- 
component dilute Bose gas [s^ |34|-|36]. Similar to the 
one-component case, we will use the RG equations 
derived at the critical point (zero particle density), to 
describe the gas of a finite but small density, and to stop 
the RG flow at some scale I — I* where the system is 
effectively no more dilute, i.e., where the "running" total 
density p{l) = ptote'*' becomes comparable with Aq. This 
condition determines the stopping scale I* as 



dl* 



tot, 



(10) 



where Cd=i = and Cd=2 = '^^^ identified by 
comparing to the known results for the one-component 
case 1361. 



III. TWO-COMPONENT BOSE MIXTURE 

Let us apply the RG approach outlined in the pre- 
vious section to the simplest case of two non-convertible 
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FIG. 1: The renormalized interaction vertex for the multi- 
component Bose gas model ^ in the zero density limit, see 
Eq. @. 



species (i.e., the particle numbers are separtely conserved 
for each species). The quartic interaction in ([1]) then 
takes the form 

C/(V'l,V'2) = Kl|^l|^+U22|V^2h/2 + Ml2|V'inV^2|', 

(11) 

SO the bare coupUngs are gn^n = wn, 522,22 = ""22, and 
512,12 = W12/2. If we denote the entries of the T matrix 
by uii, etc., the RG equations ([H]) will simplify to 



dx/dl = (2 — d)x — x'^ 



(12) 



with X e {uap}. We focus on the stability of the mixture. 
The potential (|lip is a quadratic form of the densities, 
so the necessary condition of the local stability is posi- 
tiveness of its discriminant: 



unil)u22{l) - u?2(0 > 



(13) 



which must be satisfied by the renormalized interactions 
UafsH) at all scales up to the RG stopping scale I = I*; the 
intra-species couplings Uaa must also be assumed positive 
to ensure stability. Violation of the condition (fT3t leads 
to phase separation (demixing transition) for repulsive 
interactions between species {ui2 > 0) and to a collapse 
in case of attraction {ui2 < 0). At the level of bare 
couphng {I = 0), Eq. (fT3)l corresponds to the stability 
condition obtained in the mean- field analysis [l3, El • 



In the case of one spatial dimension the solution of RG 
equations p2)l reads 



sf„MO = {i + (i/"S 



(14) 



where u^^^ denotes the corresponding bare (l — 0) value. 
A peculiarity of the ID case is the tendency of the RG 
flow to establish an enhanced symmetry if the interac- 
tion between the species is repulsive {ui2 > 0), even for 
inequivalent species. Indeed, one can see that at the fixed 
point = 1 the potential (jlip tends to 
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so the original (7(1) x (7(1) symmetry is enhanced to 
(7(1) X SU{2) at the fixed point. 

Let us first consider the repulsive case. Assume that 
at the microscopic scale {I = 0) the mixture is stable, i.e., 
< U11M22. Then the stability condition ([T3| breaks 
down starting from a certain scale I = Ic, determined by 
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The instability develops if > Ic, which translates into 
the requirement for the density to be low enough, 



Ptot < CiAoe 



(16) 



Note that Ic is not replacing the RG stopping scale I*, 
but is a new scale that naturally arises in the prob- 
lem. Here u^^'^^ denotes the mean-field result for the 
instability point [13, [lii, but the actual instability, as 
it is easy to see, generally occurs already at smaller 
values of ui2. Indeed, Eq. (|15p has a real solution Ic 
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< M12 < 
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and one has Ic 



00(0) as 



U12 — ^ (^12 )• 11 1'^^ ultraviolet cutoff Aq is much 
larger than all the other energy scales, the critical den- 
sity pc^ below which the instability develops for a given 
M12, becomes cutoff-independent and is given by 



Pc 



CiUi2{uiimi + U22'm2){ui2 — u 



(Id)-, 
12 / 



,('"/)l2 



(17) 



The above formula can be alternatively viewed as de- 
termining the critical value of the interaction u^g^ for a 
given density ptot ~ Pc- One can see that the actual 
instability point u^^2 ~^ '"12'''' if Ptot 0, and tends to 
the mean-field answer u^j' "i™^'* when the density 
becomes large. The resulting phase diagram for the one- 
dimensional case is schematically shown in Fig. [51 For 
equivalent species (un = M22, rni = TO2), the points 
u'^i^'^ and Mi™^'' coincide, and there is no change in the 
stability condition compared to the mean-field case. 

The case of attractive interspecies interactions (ui2 < 
0) is much more complicated. The RG flow for Ui2(0 
exhibits a runaway singularity at the characteristic scale 
Ih ~ In -, 1 ^° , namely, ui2(0 — ^ —00 as / — It is 

|"12| ' ^ ' 

easy to see that this scale is connected to the presence of 
molecules (bound states), which exist in one dimension 
at any strength of the attraction M12. Indeed, the energy 
of the bound state for the contact interaction U{x) = 
— |mi2|(5(x) is El, = — imi2it^2, and the total energy of 
a molecule with the total momentum Q is i^moKQ) = 
Eb + Q^/2(mi -I- 7712); molecules with momenta \Q\ > 
Qb = Aqc"''' can be broken up in collisions, while those 
with IQI < Qb are stable. 

Molecules represent a new "species" , so it is necessary 
to introduce a new molecular bosonic field and in 
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FIG. 2: (Color online). The phase diagram of a two- 
component one- dimensional dilute Bose mixture. The behav- 
ior of the critical density p in the vicinity of Ui2 — u[ 
linear, cf. Eq. (fT7|) . 
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addition to atom-molecule and molecule-molecule inter- 
action one has to include three-particle amplitudes of the 
type \E'mViV'2- Presence of those amplitudes invalidates 
the initial model ([IJ , so the attractive case requires a sep- 
arate analysis which is beyond the scope of the present 
paper. 

I would like to conclude this subsection with a re- 
mark concerning the relation between the present study 
and that by Cazalilla and Ho [23|. On the basis of the 
bosonization analysis, they have pointed out that the sta- 
bility condition is of the mean-field type for large densi- 
ties (quasicondensate regime) but changes at low densi- 
ties. Their instability criterion for the low density case 
(which corresponds to the Tonks limit rriaUaa/ Pa 3> 1, 
here pa are the densities of individual species) reads 



|m12| > 7I'^(piP2/?TllTO2)^/^, 



(18) 



and thus is different from that derived in the present 
work. However, one can argue that the condition ^TE\\ 
should be rather taken as the applicability limit of the 
bosonization approach itself. Indeed, the perturbative 
approach of [2J| can be expected to break down when the 
inter-species interaction ui2 becomes comparable with 
the characteristic bandwidths Va = i^Pa/ma of the in- 
dividual species. One can also disprove the criterion (fT8)) 
by considering the limit of hardcore bosons un, U22 ~> 00 
on a lattice: in this limit our model is formally equivalent 
to that of two fermionic species with an on-site interac- 
tion, i.e., to the one-dimensional Hubbard model. This 
model is exactly solvable and from its phase diagram it 
is known [37j that there is no instability for arbitrary 
small densities. At the same time, the lower boundary 
for the instability point ^12'^^ given by (fT5)) . diverges as 
Uaa — > 00, SO the absence of instability in the hardcore 
limit is correctly reproduced in the present approach. 



B. d = 2 

In the two-dimensional case the solutions of RG equa- 
tions are given by 

u^p{l) = {l + l/u^^l]~\ ua,pil)^'^u„p{l). (19) 

In contrast to the one-dimensional tendency to the 

enhanced SU (2) symmetry is only present for species of 
equal masses. 

Following a similar route as for d ~ 1, one obtains the 
scale Ic at which the stability condition breaks down: 
it is given by the positive root of the equation 



'mi — m2 



. mi -I- m2 
-f27r/ 



I' 



mim2 



1 



L uiimi 



1 



U11U22 
4 



U22m2 Ui2(mi+m,2) 



} (20) 



= 0. 



In the general case of unequal masses (mi 7^ m2), the 
solution If. exists for any U12 within the interval 



U12 < u]^^^] the scale Ic vanishes for 1*12 



*12 



< 
and 



diverges at U12 — > 0. An instability occurs if Ic < I* , i.e., 
if the total density is low enough to satisfy the inequality 



Ptot < C'2Aoe" 



-21^ 



(21) 



This means that, in contrast to d = 1, in two dimensions 
a dilute Bose mixture of two unequal species is always 
unstable below a certain density, for any value of the 
interspecies interaction U12, see Fig. [3Ka). However, for 
small U12 the critical density is exponentially small: ac- 
cording to Eqs. (fTS)) . (fT7|) . at U12 +0 the critical value 
of the density behaves as pc — C2^q exp(— ^/ui2), where 



A 



47r(mi -t- m2) 



y/mim2{y/rnl 



(22) 



At mi — > m2 the coefficient A diverges; actually, in that 
limit the critical density rapidly goes to zero in the whole 
interval < M12 < u'{r^\ where 



(2d) _ 



2^11^22 
"11 + "22 ' 



(23) 



Thus, the case of equal masses is special: Ic exists only 
for "12'*'' < "i™^'*, resembling the situation in the one- 
dimensional case, see Fig. ^h). If not only the masses, 
but also the intraspecies couplings of the two species are 
equal (uii = M22), the mean-field answer for the instabil- 
ity threshold is restored at all densities. 

The case of attractive interspecies interaction suffers 
from a similar runaway singularity as in the ID case. The 
scale lb = 7r(mi + m2)/mim2|wi2|, at which the running 
coupling "12(0 diverges, is connected to the existence of 
a bound state (molecule) with the binding energy Eh oc 
exp(— 2Z;,). Presence of molecules invalidates our original 
framework ([T]), so it cannot be directly applied to the 
attractive case. 
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FIG. 3: (Color online) The phase diagram of a two-component 
dilute Bose mixture in two dimensions: (a) the general case; 
(b) the special case of equal species masses mi = m2. The 
behavior of the critical values of 1/| ln(p/Ao)| in the vicinity 
of Ui2 = and ui2 — u{2'' linear, cf. Eqs. (|20p . H2ip . 




FIG. 4: RG flow diagram for the spinor boson model (|24p at 
zero magnetic field /i = 0, for the one-dimensional case. 

IV. F = 1 SPINOR BOSONS 

In this section I would like to illustrate the applica- 
tion of the presented method to a system with convert- 
ible species, where the total particle number of individual 
species is not conserved. Consider a gas of bosons with 
the hyperfine spin F = 1. The contact interaction be- 
tween two paticles depends on their total spin S*, so there 
are two characteristic scattering lengths oa with S* = 0, 
2, and the interaction can be written as [l3 . [l^ 

(24) 

where a,b & {0, ±1} denote the three components of 
the bosonic field, F^^ are the spin-1 matrices, and 5^ = 
^a^ofc^'' ^^'^ spin operators. The couplings are given 
by Co = (.90 + 252)73, C2 = (g2-5o)/3, with 55 cx as, and 
the mass m is the same for all three components. The 
last term in (f24| describes the quadratic Zeeman effect 
caused by an external magnetic field. This type of inter- 
action favors ferromagnetic spin correlations for C2 < 0, 
and polar (nematic) ones for C2 > [l^, [l^ . 

Although (1^^ includes just three parameters, the 9x9 
interaction matrix gab,a'b' contains generally five difi'erent 



couplings: 

9++,++ = fJ — . — = ri-^CQ + C2-h 

2.9+0, +0 = 2.9-0,-0 = / 1-^ Co + C2 

2.9+-,+- = V ^ Co - C2 + h, (25) 

500,00 = u^co, .9+_,oo = u"-^ 2c2, 

where arrows indicate the corresponding bare values. 

At zero external field {h = 0) the SU (2) symmetry 
dictates that the effective potential can contain only two 
constants cq and C2, and the RG equations can be cast 
into the form with x € {co + C2,co — 2c2}. The 
"tilded" variables are here defined as x — xm/Ao for 
d = 1 and x = xm/2Tr for d = 2. In one dimension 
{d ~ 1), there is a nontrivial stable fixed point (co — 
1, C2 — 0), which exhibits the enhanced SU{3) symmetry, 
similar to the double-species case where the fixed point 
was SU (2)-symmetric; the other two fixed points, (co = 
2/3, C2 = 1/3) and (cq = 1/3, C2 = —1/3) are unstable. 
The corresponding RG flow for d = 1 is shown in Fig. 2] 
Both for d = 1 and d = 2, there are two runaway flows: 
the one at C2 > co/2 marks the onset of the pairing state 
characterized by the formation of bound singlet pairs [3^ 
l40l |. while the other runaway flow at C2 < —cq {cq > 0) 
corresponds to the collapse instability when all atoms 
tend to bind into one giant ferromagnetic "drop" . Except 
of those two runaways, there are no other instabilities. 

In presence of the quadratic Zeeman term (h ^ 0) one 
cannot simply use the full matrix (j25p in the RG equa- 
tions, because a nonzero h leads to a relative shift in 
chemical potentials of different species, suppressing the 
|0) states for /i > and | ± 1) states for /i < 0, respec- 
tively. Therefore, if \h\ is large enough, one can neglect 
the suppressed species. For instance, for ft, > one comes 
back to a double-species problem for |+) and |— ) states, 
with a phase separation occuring at w > r, r > 0, i.e., 
ai Co + C2 > h > C2- At C2 > this phase separation 
corresponds to a transition between the nematic and fer- 
romagnetic states. For ft, < one effectively obtains a 
one-component gas, so nothing interesting happens. 



V. SUMMARY AND DISCUSSION 

The examples considered above show that a straight- 
forward generalization of the renormalization group ap- 
proach for dilute low-dimensional Bose gas (3^l3^[36| to 
the case of multicomponent systems allows one to obtain 
nontrivial results concerning the stability of Bose mix- 
tures. For inequivalent species, the stability conditions 
deviate strongly from the mean-field results; the latter 
are recovered at higher density. 

It is worth noting that the action ^ is quite general: 
except being applicable to a multicomponent gas in a 
continuum as well as in an optical lattice, it also arises in 
quantum magnetism problems dealing with the so-called 
Bose-Einstein condensation of magnons, induced by a 
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strong external magnetic field H in the vicinity of the sat- radius a±^ ^ cLap)- Thus, from Eq. ITSl) the critical values 



(™/) 

12 

3, the dependence of the interspecies scattering length 
ai2 on the magnetic field B is very well described by the 
following expression obtained in Ref. [49]: 



uration field i?^ |4l|. In a frustrated magnet the magnon fo^ the inter-species scattering length are a^^f> ~ 68.05 

dispersion may have degenerate minima at inequivalent , (m/) j.r \ a j- ^ t, c 

^ , ,,• , II ■ „ ■ and ol, ~ 79.7 mdependent of Accordmg to Rei. 

wave vectors, which gives rise to multiple species m an " 

effective model pt3 - )44i] . while Hg — H plays the role of 
the chemical potential. 

It would be interesting to test the predictions of the 
present paper numerically for a one-dimensional two- 
component Bose-Hubbard model with on-site interaction, 
described by the Hamiltonian 

n = -E^"(^l«Vi,a+h.c.) (26) 

j,a 

+ ^ ^ Uaanj^a{nj,a " 1) + £^12 ^j,l^j\2, 



ai2 



(27) 



is a bosonic 



where a = 1,2 labels the two species, bj^a 
operator at the lattice site j, and nj,a — -j a-j a 
the corresponding number operators. A simple estimate 
shows that the necessary requirements can be satisfied 
at realistic values of the model parameters. Taking, for 
example, two species with different hopping amplitudes 
ti = t, t2 = t/5 and the same intra-species repulsion 
t^ii = t^22 = U = 3t, from ([TS]) one obtains the follow- 
ing window in the inter-species coupling where the phase 
separation can take place: 0.745 < U12/U < 1. Further, 
fixing, for instance, U12 = Q.8U, we find from that, 
in order for the instability to develop, the total density 
per site n — ptot/Ao must satisfy n < 0.16, which can be 
reasonably reached in a density matrix renormalization 
group calculation [45|. 

Experimentally, the most promising candidate is the 
®''Rb-'*^K mixture, for which a very fine control of inter- 
actions has been achieved [81]. Assuming that all atoms 
are in their |1, 1) state (as in the setup of Ref. P), at low 
magnetic field B < 100 G the intra-species scattering 
lengths are approximately given by their B ^ values 
an ~ 97.4 for ^^Rb [H and 022 ^ 65.2 for ^iR [13, 
in units of the Bohr radius. The effective couplings in 
a cigar-shaped (quasi-lD) trap with the transverse con- 
fining frequency uj± are given by Uap = 2huj±aaj3 [3 
(where it is assumed that the characteristic confinement 



with atg ^ 284 a.u., Ai 37 G, Bi = 39.4 G, 
A2 = 1.2 G, B2 = 78.92 G. Therefore, the window for a 
phase separation a']^^^ < ai2 < 0,^12'^^ translates into the 
field range of 93.8 to 96.3 G. Taking ujj_ = 2tt x 10^ Hz, 
we obtain from Eq. that at the typical particle den- 
sity p = 10''cm~^ the phase separation instability will 
appear already at B = Be — 94.9 G, i.e., before reaching 
the mean-field threshold value Bmf — 96.3 G. Since re- 
cent experiments give hope to the possibility to control 
the interspecies scattering length 012 within a precision 
better than one Bohr radius Q , a detection of this effect 
must be within the experimental reach. 

A similar estimate can be done for the *^Rb-^^K mix- 
ture in a two-dimensional trap. The couplings in that 
case [2l[ are given by = {4:Trh!^uj±/mai3y^'^aai3. Us- 
ing a rough estimate Aq ^ l/?'Rb, where ORb — 248 pm 
is the atomic radius of Rb, one can obtain from Eqs. 
(PO)) . (PT|) that at the same transverse confining frequency 
uj± = 2tt X 10^ Hz and a typical density p — IQ^cm^^ the 
phase separation instability develops at the critical value 
of ai2 ~ 54.6, which with the help of Eq. P?)) translates 
into the critical field B ~ 91.4 G, well below the mean 
field value of 96.3 G (and this difference can be increased 
by lowering the density). 
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